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ABSTRACT 


Thin  report  presents  numerical  results  in  the  form  of  graphs  of  the 
power  reflection  coefficient  for  electromagnetic  signals  normally  inci¬ 
dent  upon  a  plasma  gradient.  An  electron  collision  frequency  that  is 
independent  of  position  in  the  plasma  interface  is  assumed.  A  "slide 
rule"  is  provided  to  enable  quick  determination  of  the  magnitude  of  the 
reflection  coefficient  for  "very  overdense,  exponential  plasmas"  in 
the  lower  atmosphere  (to  70  km)  for  frequencies  from  30  MHz  to  10,000 
MHz.  For  plasmas  that  only  reach  a  finite  density,  a  complete  set  of 

graphs  is  given  presenting  power-reflection  coefficients  for  normalized 

-20 

plasma  parameters.  Reflectivities  are  plotted  if  greater  than  10 

-3  +3 

Normalized  collision  frequencies  (Z)  from  10  to  10  are  covered. 

Normalized  plasma  to  RF  frequency  ratios  inside  the  plasma  cover  the 

-4  +4 

range  from  0.01  to  100  (10  S  X  ^  10  )  .  Characteristic  distances 

over  which  the  plasma  density  exponentiates  in  the  gradient  region  vary 
-3  +3 

from  10  wavelengths  to  10  wavelengths.  Simple  aids  to  the  reader 
are  provided  so  that  he  may  easily  determine  his  parameters  for  use 
with  the  reflectivity  curves,  should  his  plasma  be  of  atmospheric 
constituents . 
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1  INTRODUCTION 


This  report  presents  numerical  results  for  the  power-reflection 
coefficient  for  electromagnetic  waver-  normally  incident  upon  an  inter¬ 
face  between  vacuum  and  a  plasma.  Specif ically,  numerical  results  are 
presented  for  situations  in  which  the  plasma  density  varies  gradually 
over  a  wavelength  of  the  electromagnetic  w»ve .  Collision*  between  the 
plasma  electrons  and  a  neutral  background  gas  are  included  . 

Two  types  of  plasma  density  profiles  are  considered.  In  the  first, 
the  plasma  density  is  allowed  to  exponentiate  to  infinite  density.  For 
many  applications,  this  approximation  is  adequate  because  in  practice 
wave  attenuation  prevents  the  incident  signal  from  actually  interacting 
significantly  with  the  very  high  electron  densities  existing  at  great 
penetration  depths.  Therefore  the  details  of  electron  density  very  deep 
within  the  medium  do  not  matter.  The  solution  for  the  power  reflectivity 
in  this  case  is  a  very  simple  analytic  expression  that  admits  to  simple 
evaluation.  The  power-reflection  coefficient  in  a  gradient  of  this  sort 
is  presented  in  the  form  of  a  "slide  rule"  to  allow  the  reader  to  choose 
his  own  exponentiating  distance.  The  reflection  coefficient  is  presented 
as  a  function  of  altitude  in  the  earth's  atmosphere  for  frequencies  from 
30  MHz  to  10,000  MHz. 

The  second  density  profile  that  is  considered  is  one  form  of  the 
"Epstein"  profile — a  profile  that  permits,  in  principle,  a  very  wide 
latitude  in  choice  of  plasma  layer  shapes.  In  practice,  a  numerical 
solution  requires  evaluation  of  rather  unfamiliar  analytic  functions. 

For  our  example,  the  function  is  the  gamma  function  of  complex  argument. 
The  plasma  profile  considered  here  varies  monotonically  from  vacuum  to  a 
finite  electron  density.  The  electron  collision  frequency  is  assumed 
to  be  constant  throughout  the  plasma  (electron-neutral  collisions  are 
assumed  to  dominate)  .  For  very  high  electron  densities  deep  within  the 
plasma  region,  this  particular  profile  appears  to  be  nearly  the  same  as 
the  exponential  profile  considered  in  the  first  part. 
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The  determination  of  the  absolute  magnitude  of  gamma  function  of 
complex  arguments  used  in  evaluating  the  reflectivity  of  the  Epstein 
layer  requires  evaluation  of  a  continued  product.  A  computer  test  has 
revealed  that  for  many  cases  of  interest  in  this  work,  10,000  terms  of 
the  continued  product  have  not  converged  the  reflectivity  to  within  20 
percent  of  the  proper  result.  We  have  therefore  developed  a  technique 
of  analytic  approximation  to  the  continued  product  that  allows  very 
rapid  convergence  to  one-percent  accuracy.  Details  are  given  in  the 
Appendix.  The  technique  may  be  of  interest  to  others  faced  with  evalua¬ 
tion  of  similar  analytical  functions.  One  form  of  our  analytic  technique 

enables  slide-rule  calculation  of  power  reflectivity  to  an  accuracy  of 

-20 

10  percent  should  the  reflectivity  be  greater  than  10  .  Therefore  the 

reader  doe''  have  at  his  disposal  a  rather  simple,  analytical  technique 
that  he  may  apply  for  situations  of  particular  interest  not  covered  by 
our  normalized  graphs. 

Section  2  of  the  report  presents  background  information  concerning 
plasma  properties.  Section  2  also  provides  disclaimers,  reservations, 
admonitions,  and  qualifications  reflecting  the  author's  own  practical 
point  of  view  regarding  reflectivity  calculations.  Section  3  describes 
the  reflectivity  from  an  exponential  plasma  gradient  to  infinite  elec¬ 
tron  density.  Section  4  presents  the  results  of  calculations  for 
gradients  to  finite  electron  density.  The  mathematical  developments 
are  contained  in  the  Appendix. 

This  report  will  show  that  coherent  reflection  from  an  interface 
between  a  plasma  and  a  vacuum  can  be  extremely  reduced  if  the  interface 
is  gradual.  The  coherent  reflection  can  become  very  small  for  some 
plasmas  that  actually  have  very  high  electron  densities.  In  this  cir¬ 
cumstance,  it  seems  likely  that  energy  can  be  returned  to  the  s  .gnal 
source  by  any  one  of  many  other  kinds  of  reflection  mechanisms.  In 

order  to  decide  if  this  may  be  the  case,  the  author  suggests  that  if  the 
* 

Fresnel  term  for  scattering  from  the  plasma  is  large  (say,  greater 

* 

We  use  this  terminology  to  mean  the  reflectivity  that  would  result  if 
the  plasma  gradient  wore  very  short  compared  with  the  wavelength  of 
the  electromagnetic  wave.  Section  4  and  the  Appendix  further  clarify 
this  term. 
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than  10  *"  to  10  ),  yet  appllcatloi.  of  the  gradient  degrades  this  to 

-10 

less  than,  say,  10  ,  then  some  other  scattering  mechanism  Is  liable 

to  be  more  important  than  the  coherent  reflection  from  the  plasma 
gradient.  An  example  would  be  scattering  from  turbulence  or  plasma 
waves  that  produce  density  fluctuations  within  the  plasm a  gradient. 

This  scattering  is  best  calculated  using  Born  approximation  and  a  sto¬ 
chastic  description  of  the  dielectric  fluctuations.  Furthermore,  under 
these  same  conditions  of  large  Fresnel  term,  incoherent  backscatter  from 
free  electrons  will  begin  to  be  as  intense  as  coherent  effects  due  to 

the  mean  plasma  structure  at  power  reflectivity  values  below  the  order 
-20 

of  10  For  this  reason,  computer  plotted  data  presented  here  do  not 

-20 

extend  below  10 
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2  BACKGROUND 


The  complex  dielectric  content  of  a  plasma  is  described  approxi¬ 
mately  by  the  Appleton-Hartree  formula.1*  Under  the  assumption  that 
magnetic  field  effects  are  small,  this  dielectric  constant  is  given  by 


C  =  1 


(1) 


where 

X  =  (f  /f)2 
P 

f  =  Plasma  frequency  in  Hertz 

f  =  Frequency  of  the  electromagnetic  wave  in  Hertz 

Z  =  V/UD  =  f  ,/f 
c 

v  =  Collision  frequency  in  radians/second 
u)  =  2nf 

f  =  Collision  frequency  in  Hertz. 

The  plasma  frequency  is  related  to  electron  density  through 

f  =  8.97  •  103/IT  Hertz  (2) 

P  e 

where  ne  =  electron  density  in  number  per  cubic  centimeter. 

In  the  earth's  atmosphere,  Shkarofsky3  has  shown  that  the  quantities 

Z  and  X  used  in  Eq .  (1)  are  not  quite  as  we  have  defined  them.  The 

energy  dependence  of  the  electron  collision  cross  soction  with  nitrogen 

is  such  that  one  must  use  effective  values  for  X  and  Z  that  are  related 

to  our  definitions  in  a  complex  way.  Figures  1  and  2  present  curves 

for  relating  our  definition  of  X  and  Z  to  X  and  Z  to  be  used  in 

eff  eff 

calculations.  For  the  data  of  Sec.  3,  we  have  not  made  this  correction, 

* 

References  are  listed  at  vne  end  of  this  report. 
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though  the  effect,  in  a  practical  sense,  is  very  minor.  The  reader  may 
do  as  he  wishes  in  using  the  data  of  Sec.  4. 

To  aid  the  reader  we  present  in  Figs.  3  and  4  our  interpretation  of 
the  proper  electron-neutral  collision  frequency  as  a  function  of  alti¬ 
tude  in  the  earth's  atmosphere.  In  Fig.  3  the  collision  frequency  is 
presented  in  radian  measure.  In  Fig.  4  the  collision  frequency  is  pre¬ 
sented  in  megahertz  in  order  that  the  reader  may  more  readily  compute  Z 
for  his  radar  frequency.  Equation  (2)  is  plotted  in  Fig.  5;  this  curve 
gives  the  plasma  frequency  in  megahertz  versus  electron  density,  which 
is  needed  to  compute  X . 

The  calculations  to  be  presented  assume  that  the  electron  collision 
frequency  is  constant  throughout  the  plasma  gradient.  This  situation 
pertains  to  a  plasma  in  which  the  electron-neutral  collision  frequency 
is  larger  than  other  types  of  collision  frequencies.  The  electron-ion 
collision  frequency  is  proportional  to  electron  density.  Therefore,  as 
the  electron  density  increases,  electron  collisions  with  ions  may  be¬ 
come  as  large  as  the  electron-neutral  collision  frequency.  In  the  plasma 
density  gradient,  then,  electron-ion  collisions  may  begin  to  dominate  the 
local  propagation  parameters.  Our  analytical  results  could  be  in  error 
under  these  circumstances. 

Figure  6  presents  the  electron-ion  collision  frequency  versus  elec¬ 
tron  density  for  three  temperatures  of  air.  Along  the  right-hand  ordinate 
we  have  indicated  altitudes  from  Fig.  3  at  which  the  electron-neutral 
collision  frequency  attains  the  same  value  as  the  ordinate  value.  The 
reader  must  use  care  in  employing  the  results  to  be  presented  later, 
particularly  for  high-frequency  radars.  That  is,  if  the  electron  density 
in  his  plasma  profile — for  his  air  temperature — produces  an  electron 
collision  frequency  that  equals  or  exceeds  at  some  depth  the  electron 
neutral  collision  frequency,  then  he  must  attempt  to  determine  whether 
penutration  of  his  waves  to  that  depth  significantly  affects  the  re¬ 
flectivity  to  be  expected.  The  reflectivity  may  be  increased  or  de¬ 
creased  by  this  effect,  according  to  details  of  the  situation.  We  have 
not  studied  reflections  from  electron-density  profiles  in  which  collision 
frequency  varies  with  position. 
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ELECTRON -NEUTRAL  COLLISION  FREQUENCY - rod/sec 


ELECTRON -NEUTRAL  COLLISION  FREQUENCY 


PLASMA  FREQUENCY - MHz 


ELECTRON  COLLISION  FREQUENCY - tob/wc 


The  data  to  be  presented  were  calculated  for  plasma  gradients  that 
are  expressed  in  term:  of  simple,  analytic  functions.  It  is  hardly  ex¬ 
pected  that  real  plasma  gradients  will  take  on  the  forms  that  we  use. 

We  suggest  therefore  that  our  data  be  used  primarily  as  a  guide  as  to 
whether  a  coherent  plasma  interface  reflection  or  some  other  reflection 
or  scattering  mechanism  dominates  a  given  measurement.  We  also  suggest 
that  if  the  plasma  profile  expected  is  very  poorly  approximated  by  our 
analytic  expressions  but  the  plasma  frequency  inside  considerably  ex¬ 
ceeds  the  radar  frequency,  then  the  reader  should  match  our  analytic 
form  to  his  plasma  gradient  at  a  position  in  the  gradient  where  the 
radar  wavelength  in  the  plasma  is  longest  compared  with  the  plasma 
gradient.  We  can  show  that  the  radar  wavelength  in  a  plasma  is 


where 

\  =  Radar  wavelength  in  plasma 

=  Radar  wavelength  in  free  space. 

The  wavelength  takes  on  its  greatest  value  at  X  =  2,  if  Z  is  finite. 
If  Z  is  zero,  then  the  power  reflectivity  is  1.0  regardless  of  plasma 
profile;  here,  profile  details  would  be  of  concern  only  if  the  absolute 
phase  of  the  return  signal  were  of  interest.  The  limit  of  \  with  Z  =  0 
is  never  met  in  practice,  however,  though  some  environments  do  approach 
this  situation.  Therefore,  in  a  situation  where  the  reader  encounters 
a  physical  plasma  profile  very  different  from  our  analytical  forms,  he 
might  match  our  analytic  gradient  at  the  position  where  X  =  2  (f^  =  /2  f) . 
On  the  other  hand,  our  intuition  suggests  that  penetration  of  the  wave 
to  the  X  =  2  depth  may  not  be  important,  so  that  the  gradient  matching 
should  bo  done  at  the  X  =  1 .0  depth,  particularly  for  very  small  values 
of  Z.  The  choice  for  matching  is  left  to  the  reader. 

The  main  effect  of  electrons  in  the  region  where  X  <  1.0  is  probably 
wave  attenuation.  Therefore,  should  the  electron  density  protile  in  this 


region  be  very  different  from  our  analytic  forms,  the  reader  may  still 
wish  to  locally  fit  his  physical  gradient  to  our  analytic  gradient  in 
the  region  from  X  =  1  to  X  =  2 .  He  would  then  compute  by  WKB  approxi¬ 
mation  the  difference  in  absorption  in  the  region  where  the  profiles  are 
very  different,  in  order  to  obtain  a  more  appropriate  reflectivity. 

Perhaps  the  most  useful  application  of  the  data  presented  here  would 
be  to  check  results  of  computer  programs  that  are  designed  to  determine 
plasma  reflectivities  using  a  layer-by-layer  numerical  mock-up  of  the 
physical  environment.  One  extreme  danger  of  the  mock-up  procedure  is 
that  the  amount  of  computer  timb  required  increases  with  the  number  of 
layers  considered.  Since  the  total  number  of  layers  that  can  be  handled 
may  be  limited  by  computer  capacity,  the  layer  thicknesses  tend  to  be 
picked  too  large.  As  a  result,  at  each  interface  between  layers  (or  be¬ 
tween  shells  for  a  spherical  geometry),  the  reflectivity  may  be  larger 
than  the  total  reflectivity  of  the  medium.  The  person  performing  this 
kind  of  calculation  then  expects  that,  since  he  will  be  preserving  phase, 
reflectivities  from  all  the  layers  will  tend  to  phase-cancel  properly 
to  leave  the  correct  residual  reflectivity.  The  difficulty  with  this 
procedure  is  that  the  reflectivity  of  one  interface  is  likely  to  domi¬ 
nate,  while  all  other  interface  reflectivities  simply  produce  a  sort  of 
random  walk  away  from  this  dominant  interface  value. 

Our  analytical  results  can  be  used  to  check  such  computer  programs. 
The  procedure  would  be  to  mock-up  our  analytic  form,  layer  by  layer, 
compute  reflectivities,  and  then  compare  results.  Good  agreement  would 
help  develop  confidence  in  the  layer-by-layer  mock-up  procedure.  Then 
the  mock-up  procedure  could  be  applied  to  both  plasma  gradients  and 
collision  frequency  variations  that  cannot  be  handled  analytically. 
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3  REFLECTION  FROM  AN  EXPONENTIAL  PLASMA  GRADIENT 


The  mathematical  derivations  of  reflectivities  used  in  this  report 
were  taken  from  Budden  .3  Section  17.2  of  Ref.  3  derives  the  reflectivity 
to  be  expected  in  a  plasma  gradient  in  which  the  electron  density  varies 
exponentially.  The  derivation  assumes  uniform  collision  frequency.  We 
def ine 

z/z 

n  (z)  =  n  e  °  (4) 

e  o 


whore 

ne(z)  =  Electron  density  at  position  z 

=  Electron  density  at  position  z  =  0 

z  =  Position  in  gradient  ( -°°<  z  <  +  ®) 

z^  =  Length  over  which  the  electron  density  changes  by  a 
factor  of  o.  This  value  would  be  matched  to  the 
reader's  problem. 


In  terms  used  in  the  previous  section,  Eq .  (4)  becomes 

z/z 

X  «  c  ° 


(5) 


For  normal  incidence,  Ref.  3  gives,  for  the  power  reflection  coefficient 
[Eq.  (17.17)  of  Ref.  3], 


exp  j- 


81  Tz 


o  -1 
—  ton  Z 


(6) 


Figures  7  and  8  are  to  be  used  to  evaluate  Eq  .  (6).  Figure  7  is  a 
transparency  to  be  used  as  an  overlay  on  Fig.  8.  The  ordinant  of  Fig. 

8  gives  the  power  that  10  must  be  raised  (always  negative)  in  order  to 
obtain  the  reflectivity. 
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FIG.  7  OVERLAY  FOR  USE  WITH  FIG. 
FROM  EXPONENTIAL  PLASMA 


8  FOR  CALCULATION  OF  REFLECTIVITY 
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REFLECTIVITY  EXPONENT 


EXPONENTIATING  DISTANCE 


An  acetate  version  of  Fig.  7  (in  pocket  under  back  cover)  is  over¬ 
laid  on  Fig.  8  so  that  the  arrow  labeled  z  is  adjacent  to  the  gradient 

o 

distance  of  interest.  For  an  exponential  gradient  of  that  value,  the 
reader  obtains  the  reflectivity  exponent  for  any  altitude  or  frequency 
of  concern.  As  an  example,  let  the  exponentiating  distance  be  one  meter; 
then  at  an  altitude  of  10  km,  the  power-reflection  coefficient  for  a 

-53 

1000-MHz  radar  wave  from  a  very  overdense  plasma  is,  approximately,  10 

clearly  a  reflectivity  of  little  concern.  On  the  other  hand,  a  100-MHz 

—5  8 

radar  wave  will  experience  a  reflectivity  value  of  10  ,  which  may  be 

of  considerable  concern. 

At  an  altitude  on  the  order  of  40  km,  at  a  plasma  frequency  of 
10,000  MHz  and  a  plasma  temperature  of  1000>'K,  the  electron-neutral 
collision  frequency  and  the  olectron-ion  collision  frequency  are  com¬ 
parable.  At  higher  altitudes  than  this  the  electron-ion  collision 
frequency  is  greater  than  the  electron-  leutral  collision  frequency. 
Therefore,  for  10,000  MHz,  a  plasma  temperature  of  1000’K,  and  altitudes 
above  40  km,  the  reflection  coefficient  is  not  given  by  Eq .  (6).  Wo 
have  shown  a  shaded  region  that  we  believe  probably  contains  the  proper 
reflectivity  value;  at  this  writing  we  are  not  prepared  to  provide  a 
precise  value. 

Similar  shaded  regions  are  given  for  6000  and  3000  MHz,  and  1000 JK 
temperatures.  For  higher  temperatures,  the  'lectron-ion  collision  fre¬ 
quency  is  lower,  so  that  greater  reflectivities  will  be  obtained.  We 
have  chosen  1jOO°K  for  this  discussion  because  air  temperatures  would 
generally  have  to  be  at  least  this  high  in  order  to  maintain  electron 
densities  by  radiation  (or  thermally)  that  are  overdense  to  frequencies 
of  3000  MHz  and  above. 


18 


4  REFLECTION  FROM  A  GRADUAL  GRADIENT  TO  FINITE  ELECTRON  DENSITY 


In  this  section  graphs  are  presented  that  give  the  power  reflection 
coefficient  for  electromagnetic  signals  normally  Incident  upon  a  plasma 
cl  finite  electron  density.  The  electron  density  In  these  cases  varies 
from  a  value  of  zero  per  cubic  centimeter  very  far  from  the  plasma  to  a 
value  of  ciQ  deep  Inside  the  plasma.  The  analytical  form  for  the  varia¬ 
tion  of  electron  density  Is 


n  (z) 
e 


1 


n 

o 


+  e 


-z/a 


(7) 


The  parameter  a  characterizes  the  steepness  of  the  plasma  gradient. 
Vacuum  exists  at  z  =  -  ®;  plasma  with  electron  density  n^  exists  at 
z  =  -*•  ®.  For  very  large  negative  values  of  position  z, 


n  (z) 
e 


+z/a 

n  e 
o 


(8) 


This  analytical  form  is  like  that  of  Sec.  3.  Should  nQ  be  large  enough 

so  that  the  wave  does  not  penetrate  to  a  depth  where  Eq .  (8)  Is  a  poor 

approximation  to  Eq .  (7),  then  results  of  this  section  will  reduce  to 

those  of  Sec.  3  with  o  =  z  of  Sec.  3. 

o 

The  plasma  dielectric  constant  in  the  plasma  gradient  is  given  by 


C(z) 


1 


X 

1  -  1Z 


1  +  e 


-z/a 


(9) 


The  collision  frequency  for  this  form  is  constant  as  a  function  of  posi¬ 
tion  in  the  plasma  gradient.  This  condition  approximates  plasmas  in  the 
lower  atmosphere.  The  power-reflection  coefficient  for  this  analytical 
gradient  form,  given  in  Ref.  3,  is 


1  - 

2 

r Cl  +  ikc(/i^  +  1)  ] 

1  +  /«! 

r  [1  +  ikr-(/G^  -  1)  ] 
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whore  Cg  is  the  dielectric  constant  from  Eq .  (9)  at  z  equal  to  plus 
infinity.  The  parameter  k  is  equal  to  2tt/X.  The  first  factor  on  the 
right-hand  side  of  Eq .  (10)  is  the  familiar  Fresnel  equation  for  power 
reflection  from  an  abrupt  interface  between  vacuum  and  plasma.  For  this 
leason  we  refer  to  this  factor  as  the  Fresnel  term  in  this  report.  The 
second  factor,  involving  the  gamma  function  ratio,  introduces  the  effect 
of  the  change  in  electron  density  occurring  over  a  finite  distance.  It 
is  clear  that  if  o  goes  to  zero,  meaning  an  infinitely  sharp  transition, 
then  Eq .  (10)  is  left  with  only  the  Fresnel  term,  as  is  necessary.  The 
Appendix  indicates  in  detail  how  Eq .  (10)  was  evaluated  and  the  accuracy 
that  we  expect  in  our  plots. 


The  power-reflection  coefficient  is  plotted  as  a  function  of  Z, 

f  /f,  and  o/X  in  Figs.  9  and  10.  In  Fig.  9(a)  through  9(m)  the  logarithm 

of  the  power-reflection  coefficient  is  plotted  versus  the  logarithm  of 

the  ratio  o/X.  Each  graph  page  is  for  a  fixed  value  of  Z  for  values  of 
-3  +3 

Z  from  10  to  10  .  The  power-reflection  coefficient  was  computed  and 

plotted  versus  o/X  for  a  fixed  f  /f  ratio,  where  f  is  the  plasma  fre- 

quoncy  and  f  is  the  radar  frequencv.  The  choices  of  values  for  this 

ratio  that  were  plotted  were  made  for  each  graph  so  as  to  minimize  gaps 

between  adjacent  curves.  Thus,  different  choices  for  the  values  of  f  /f 

*  P 

will  appear  on  various  graphs.  The  power-ref lect ion  coefficient  was 

2 

computed  for  Z  and  a  fixed  value  of  X,  which  is  (f  /f)  .  The  confuta¬ 
tion  was  performed  for  values  of  a/X  that  were  stepped  in  the  ratio  of 
(10)1  in  this  way,  ten  points  were  plotted  at  even  intervals  for 

each  power-of-10  change  along  the  logarithmic  abscissa. 


The  general  behavior  of  all  of  these  curves  of  power  reflectivity 

versus  o/X,  as  presented  in  Fig.  9,  is  similar.  The  power  reflectivity 

is  constant  at  the  Fresnel  value  as  a  function  of  o/X  for  small  values 

of  j/X.  Thon  the  reflectivity  proceeds  to  decrease  extremely  rapidly 

as  the  2 Mi  ratio  continues  to  increase.  For  fairly  low  values  of  Z  and 

large  values  of  o/X,  we  note  that  the  power  reflectivity  changes  from 

almost  no  reflectivity  for  f  /f  “*  0.9,  to  a  rather  substantial  reflec- 

P 

tivitv  for  f  /i  1.1.  This  extreme  change  in  reflectivity  with  f  /f 
P  P 
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FIG.  10<m) 


for  large  values  of  o/\  aids  in  producing  relatively  unambiguous  iono- 
grams  obtained  with  Ionospheric  sounding  equlpnent.  If,  for  Instance, 
the  lonoaphere  produced  a  very  sharp  Interface  (aA  «  1.0),  then  the 
power  reflectivity  would  begin  to  decrease  slowly  with  an  Increase  In 
frequency  above  the  layer  critical  frequency  so  that  the  observed,  abrupt 
termination  of  reflection  at  F^fg,  *or  exanPle>  would  not  be  obtained. 

The  second  group  of  figures,  Fig.  10(a)  through  10(m),  presents  the 

power  reflectivity  plotted  versui  Z  for  a  variety  of  values  of  oA.  The 

power  reflectivity  displayed  in  this  manner  demonstrates  some  interesting 

behavior.  Most  startling  is  that  demonstrated  in  Fig.  10(g)  for  the 

situation  with  oA  =  1.0.  For  a  value  of  f  /f  =  0.9,  note  that  the 

P 

power  reflectivity  Increase!  very  gradually  as  Z  incre  jes  rather  than 
decreasing  as  one  would  normally  expect!  The  author  hopes  that  this  be¬ 
havior  results  from  a  physical  effect  rather  than  from  an  as  yet  un¬ 
realized  difficulty  with  his  analytics  and  approximations.  We  note  that 
a  similar  behavior  occurs  to  a  lesser  extent  in  Fig.  10(f)  for  a  aA 
ratio  of  3  X  10  * . 

Assuming  that,  indeed,  the  results  are  correct  and  the  reflectivity 
does  increase  slightly  with  an  increase  in  Z,  we  then  believe  that  the 
slight  increase  results  from  changes  in  effective  radar  wavelength. 

Thus,  at  the  proper  collision  frequency,  the  radar  wavelength  is  "in 
resonance"  with  the  plasma  gradient.  That  is,  the  currents  induced  by 
the  incident  signal  are  phased  just  right  so  that  their  reradiation  from 
the  regions  of  steepest  curvature  in  dielectric  constant  reinforce  in 
the  backscatter  direction.  At  higher  collision  frequencies,  the  wave¬ 
length  in  the  plasma  is  too  short.  At  lower  collision  frequencies,  the 
wavelength  in  plasma  is  too  long . 

For  more  shallow  gradients  (larger  oA  ratios),  f  /f  ratios  would 

P 

have  to  be  larger  than  0.9  in  order  for  the  effective  wavolength  in  the 
plasma  to  be  long  enough  to  be  "synchronous"  with  the  spacing  between 
regions  of  maximum  dielectric  curvature.  In  these  less  iteep  gradients, 
though  the  signal  wavelengths  may  be  nearly  in  resonance  for  f  /f  close 
to  1.0,  absorption  effects  could  dominate  over  resonance.  For  oA  ratios 


73 


leas  than  1.0,  the  spacing  of  regions  of  Baxleys  curvature  may  be  less 
than  the  free-space  radar  wavelength,  ao  that  no  resonance  can  occur 
(wavelenf  ths  only  Increase  In  plasaaa) . 

Our  Intuition  tells  ua  that  the  observed  lncreaoe  In  reflectivity 
will  be  mathematically  observable  only  for  plasmas  with  o/X  values  near 
1.0  and  fp/f  ratios  nearly  equal  to,  but  leas  than,  1.0.  We  suspect 
that  the  reflectivity  values  in  Fig.  10(g),  where  this  "resonance" 
appears,  are  so  very  low  that  the  coherent  reflection  from  the  plasma 
Interface  phenoaenon  that  we  are  studying  in  this  report  la  unimportant. 
Therefore  the  resonance  phenoaenon  is  of  no  practical  lnportance . 


The  second  interesting  behavior  exhibited  in  Figs.  10(e)  and  10(f) 
is  the  "plateau"  in  reflectivity  at  large  values  of  Z  and  at  values  of 
f  /f  that  are  very  large  compared  with  1.0.  We  explain  this  phenomenon 
as  follows.  For  fixed  o/\  value,  very  snail  Z,  and  very  large  f^/f 
(overdense,  collisionless  plasma),  the  reflectivity  takes  on  the  Fresnel 
value — which  is  near  1.0.  The  radar-wave  r<  flection  occurs  primarily 
in  the  region  where  f  /f  =  1.0.  As  the  collision  frequency  is  increased, 
tho  reflectivity  value  begins  to  decrease  as  the  coll  is  ion -produced  real 
part  of  the  propagation  para«§ier  allot*  t  energy  to  propagate  into  the 


overdense  plasma  region. 


For  very  large  f^/f  values 


> 


the  reflectivity 


value  will  decrease  as  Z  (the  electron  collision  frequency)  increases 


with  values  that  are  the  sane  as  those  computed  in  Sec.  3. 


As  Z  further  increases,  the  plasma  continues  to  appoar  as  an  expo¬ 
nential  plasma  at;  discussed  in  Sec.  3  with  the  arctangent  Z  term  of  Eq . 
(6)  saturated  near  m/2;  further  increases  in  Z  cause  no  observable  change 
in  arctangont  Z.  Thus,  the  reflectivity  value  level;:  off  is  a  function 
of  increasing  Z.  In  the  plasma,  as  Z  is  increased,  the  effective  region 
where  reflection!  occur  penetrates  deeper  and  deeper  into  the  medium 
even  though  the  power  reflectivity  value  remains  constant.  Eventually, 
for  u  largo  onough  value  of  Z,  the  wave  penetratoi  deeply  enough  into 
the  plasma  so  that  the  Eq .  (8)  approximation  to  Eq .  (7)  is  not  valid 
whnro  reflection  effectively  takes  place.  That  is,  for  the  exponential 
plasma  of  Sec.  3  the  region  of  oflcctive  reflection  occurs  at  a  dopth 


where  the  electron  density  exceeds  that  which  our  particular  finite- 
density  plasma  achieves.  Therefore,  our  finite-density  plasma  becomes 
ineffective  as  a  reflector,  so  that  further  increases  in  Z  are  accom¬ 
panied  by  a  decrease  in  power-reflection  coefficient.  This  intuitive 
explanation  for  the  observed  plateau  values  for  reflectivity  is  supported 
by  the  fact  that  reflectivity  values  for  the  exponential  plasma  of  Sec. 

3  agrees  with  the  plateau  values.  Figure  11  presents  plots  of  reflec¬ 
tivity  for  the  profile  forms  of  Eqs .  (7)  and  (8). 
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FIG.  11  COMPARISON  OF  POWER  REFLECTIVITY  vs.  Z  FOR  A  GRADUAL  GRADIENT 
AND  AN  EXPONENTIAL  PROFILE 


5  CONCLUSIONS 


In  this  report  we  have  presented  data  to  allow  the  reader  to  esti¬ 
mate  the  power  reflectivity  of  an  electromagnetic  signal  normally  inci¬ 
dent  upon  a  plasma  gradient.  The  electron  collision  frequency  is  assumed 
constant  throughout  the  plasma.  We  have  provided  graphs  that  allow  for 
simple  determination  of  plasma  characteristics  in  the  atmosphere.  We 
have  presented  in  Sec.  3  a  "slide  rule"  to  allow  easy  estimation  of  the 
power  reflectivity  from  an  exponential  plasma  (electron  density  goes  to 
infinity)  in  the  earth's  atmosphere.  In  Sec.  4  we  have  presented  data 
for  a  specific  analytical  form  for  the  plasma  gradient  to  a  finite, 
uniform  electron  density. 

For  reflection  from  atmospheric  plasma  gradients,  we  suggest  the 
following  procedure.  Determine  whether  or  not  the  atmospheric  plasma 
is  overdense  to  the  electromagnetic  signal  of  concern.  If  it  is,  then 
use  the  slide  rule  of  Figs.  7  and  8  of  Sec.  3  to  determine  the  power 
reflectivity — provided  the  form  of  the  gradient  is  well  matched  by  the 
exponential  form.  If  the  reflectivity  is  large  enough  to  be  of  concern, 
then  the  reader  would  be  well  advised  to  pursue  the  problem  in  some 
additional  depth  by  means  of  the  data  presented  in  Sec.  4.  The  power 
reflectivity  as  given  in  Sec.  3  will  always  be  larger.  If  ;he  reflec- 

r  -10 

tivity  is  very  small  [less  than  10  say],  then,  as  indicated  earlier 
in  this  report,  some  other  reflection  mechanism  may  be  more  important 
than  the  coherent  effect  of  the  plasma  gradient.  The  reader  would  then 
be  well  advised  to  study  other  possibilities. 

If  the  atmospheric  plasma  is  not  overdense  to  the  signal  of  con¬ 
cern,  then  the  data  of  Sec.  4  must  be  used.  Again,  the  reader  should 
be  alert  to  the  physics  of  his  environment  inasmuch  as  mechanisms  other 
than  smooth  plasma  gradients  can  reflect  electromagnetic  signals. 

In  either  case,  overdense  or  underdense  plasmas,  the  reader  can 
probably  match  well  enough  the  plasma  gradient  in  the  reflecting  region. 
Differences  in  plasma  profiles  in  the  exterior  can  be  handled  by  comparing 
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signal  absorption  by  WKB  approximation  for  the  twe  external  regions. 
Appropriate  corrpensat ion  can  then  be  made. 


i 


1 
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Append  ix 


REDUCTION  OF  ANALYTIC  FORM  FOR  REFLECTION 


FROM  A  GRADIENT  IN  PLASMA  DENSITY 


Append ix 

REDUCTION  OF  ANALYTIC  FORM  FOR  REFLECTION 
FROM  A  GRADIENT  IN  PLASMA  DENSITY 


Reference  3  (Eq.  17.112)  gives  for  the  voltage  reflectivity  between 
a  vacuum  and  a  plasma  the  following  form: 


C  -  q, 

R  =  - - - 

v  C  +  qr 


(-  2 ik 3C)  ! 
(2ik-iC) 


[iko(q, 

_ t 

[iko(qr 


+  o}:f 
■^cTTtJ 


(A —1 ) 


where 


C  =  cosine  of  incidence  angle  (C  =  1  for  normal  incidence) 
2 

q g  =  e2  complex  dielectric  constant  deep  inside  plasma 


with  X  and  Z  defined  in  Sec.  2,  and 

k  =  2tt/\ 


K  =  Radar  wavelength 

a  =  Distance  characteristic  of  change  in  plasma  density  [see 
Eq.  (A-2) ] . 

The  functions  indicated  are  factorial  functions  that  are  also  gamma 
functions.  They  have  complex  arguments. 

The  expression  (A-l)  gives  the  voltage  reflection  coefficient  (in¬ 
cluding  phase)  from  a  plasma  gradient  in  which  the  dielectric  constant 
varies  in  the  following  way: 


€(z)  =  1 


X  1 

1  -  iZ  ‘  -z/ov 

(l+o  ) 


(A-2) 


This  form  means  that  the  collision  frequency  Z  is  not  a  function 
of  position,  but  X  is .  In  fact 


n  (z)  I  n 
e  o 


(A-3) 
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For  large  negative  values  of  z  we  find  that  the  electron  density 
is  nearly  a  simple  exponential  as  was  used  in  Sec .  3.  For  very  large 
negative  values  of  z,  the  dielectric  constant  approaches  the  constant 
value  1.  At  very  large  positive  values  of  z  the  electron  density 
approaches  a  constant  value  of  nQ  and  the  dielectric  constant  similarly 
approaches  a  constant  value. 

We  shall  be  substituting  a  value  1.0  for  C  (we  ore  only  interested 
in  normal  incidence).  The  first  factor  on  the  right-hand  side  of  Eq . 

( A — 1 )  wo  call  the  Fresnel  term.  This  factor  givos  the  plasma  reflection 
coefficient  should  the  plasma  gradient  be  exceedingly  abrupt  compared 
with  the  wavelength. 

The  second  factor  on  the  right-hand  side  of  Eq .  (A-l)  has  modulus 
1.0  and  only  affects  tho  phase  of  the  returning  signal.  We  shall  be 
concerned  only  with  power  reflectivity,  so  that  this  factor  is  of  no 
concern.  However,  as  an  exercise,  for  real  o  (the  only  kind  wo  consider). 


(-  2  ik  j)  I 

2 

r(l  -  2ikj) 

2 

r<7) 

(+  2iki) : 

r ( 1  4  2 ika) 

mr 

with  z  hero  meaning  the  complex  number  (1  +  2ikj)  .  If  wo  apply  the 
ident ity 


r(z)  =  rm 


( A  —5 ) 


we  f ind 


r (f)  2  r(z)n?)  rtzin?)  _  1  n 

T(z)  T(7)?(7.)  nWruT 


(A  -6) 


An  evaluation  of  the  third  factor  on  the  right-nand  side 
straightforward.  Wo  find 


|’(iko(q2  +  '."j2 

2 

rru  +  ika(q2  I  1 )  )"]2 

[lik  (q2  -  1 )  1 t  J 

1  T  [1  ♦  iko(q2  —  1 >  1 J 

is  not 


(A-7) 
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We  let 


q2  =  M  -  iX 


and  f  ind 


1  2 
{r [1  +  kerx  +  ika(n  +  1)  ]1 

to 

tf(x  +  iy^)  ]2 

Tr [1  +  kn\  +  iko(n  -  1)  l)2 

ii 

[f(x  +  iy2)]2 

with 


x  =  1  +  kox 


y1  =  ka(u  +  1) 
y2  =  ka(n  -  1) 


(A-8) 


(A-9) 


(A-10) 

(A-ll) 

(A-12) 


Reference  4  [Eq .  (6.1.25),  p.  256]  allows  us  the  following  trans¬ 
formations  : 


or 


n! 


TT— 


n=0 


1  + 


(x  +  n)‘ 


'n — *- 


n=0 


1  + 


(x  +  n)‘ 


[ko<*  -  1)J2  | 

+  - - - 


TT  (,  ,  •  niM  -i  |i  *  [“^]  ( 


TT 

n  0 


[n  ♦  1  +  k  -x  ] 


( A-l 3) 


( A  - 1 4 ) 


We  have  found  that  the  evaluation  of  the  continued  product  of  Eq  . 
(A-ll)  may  require  tens  of  thousands  of  factors  before  the  continued 
product  converges  to  within  a  few  percent  of  the  proper  answer.  W'o 
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later  develop  a  technique  that  combines  continued  product  evaluation  and 
analytical  approximation  to  obtain  an  at  curate  representation  of  Eq  . 

(A  — 14 )  using  a  very  reasonable  amount  of  computer  time.  We  return  to 
this  shortly. 

2 

The  determination  of  |i  and  x  from  is  fairly  straightforward.  We 
find  that  by  defining 


and 


6=1 - ^—5-  (A -15) 

r  ,  „2 


V7 

«,  =  - ,  (A-16) 

1  +  Z 


and 


Kq. 


(A-17) 


(A-18) 


In  terms  of  ju  and  x,  the  first  factor  on  the  right-hand  side  of 
( A — 1 ) ,  when  absolute  squared,  becomes  the  very  familiar 


R 


F 


2  2 

Cl  -  n)  ♦  X 

x2  2 

(1  +  M)  +  X 


(A-19) 


which  is  the  power-reflection  coefficient  for  an  abrupt  interface  be¬ 
tween  the  plasma  and  a  vacuum.  This  is  known  as  the  Fresnel  equation 
and  we  call  it,  here,  the  Fresnol  term. 

We  see  from  Kq  .  (A-l)  and  (A-14)  that  should  the  parameter  j 
approach  zero  (an  abrupt  interface),  then,  since  Eq .  (A-14)  takes  the 
value  1.0,  Kq  .  (A-l)  does  reduce  to  the  expected  form. 
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Because  of  the  slow  convergence  of  Eq.  (A-14)  we  have  developed  an 
analytic  approximation  to  evaluate  the  continued  product.  In  our  compu¬ 
tations  we  have  usually  used  a  hybrid  form,  as  will  be  described.  We 
def ine 


a  =  ko(n  -  1) 

b  =  ko(n  +  1) 


(A-20) 


(A-21) 


and 


y  =  rn  +  k-ix 

m 


(A-22) 


Then  Eq .  (A-14)  becomes 


P  = 


n 

m=l 


2 

2 

y 

+  a 

m 

2 

J2 

ym 

+  b 

(A-23) 


or 


or 


P  =  exp 


!».fi  ■  ■'!  ■ 


exp 


m=l 


2  ^  |  log  [  ] 
m-1 


P  =  exp  {2 


X]  l0ge  (ym  *  *2)  l0g»  (v»  *  b2| 


(A-24) 


(A-25) 


m^l  m-1 

Translorming  the  summations  to  integrals  we  find 

C*>  00 

I  log  |y2  +  a  ’ jdm  -  I  -og  £y2  ♦  b2jdm 

Jm  1/2  Jm=l/2 


P  =  exp 


(A-26) 


The*  integrals  may  bet  ovaluatod  to  give 


P  =  oxi 


tr(a  -  b)  -  y  log 
m  ©l 


2  2 

ym  *  * 

2  v2 

y  +  b 


+  2lb  tan 


-1  ym 


-  a  tan 


-1 


(A -27) 

We  have  freedom  in  tho  way  we  may  use  this  analytical  expression 

in  that  y  ,  the  lower  limit,  can  be  chosen  at  will.  Thus,  if  we  wish 
m  7 

to  use  Eq  .  (A-27)  as  a  complete  substitute  for  the  continued  product, 
then 


ym  =  2  + 


(A-28) 


We  use  tho  1/2  in  Eq .  (A-28)  since  the  value  for  m  =  1  of  tho  con¬ 
tinued  product  is  approximated  by  the  integral  from  m  =  1 '2  to  m  =  3/2. 
Test  cases  run  on  computers  have  shown  that  Eq  .  (A-27)  may  be  in  error 
as  a  complete  representation  of  the  continued  product  of  Eq .  (A-14)  by 

as  much  as  10  percent,  but  usually,  particularly  for  total  reflect ivitie 
-20 

larger  than  10  ,  has  an  accuracy  exceeding  one  percent.  Greater 

accuracy  can  bo  obtained  by  evaluating  tho  first  m  factors  of  the  con¬ 
tinued  product,  and  then  multiplying  these  by  Eq .  (A-27)  with 


y.  *  "  4  2  *  |M* 


(A-99) 


Table  A -1  presents  a  matrix  of  values  for  power  reflectivity  in 

which  the  hybrid  form  of  m  factors  of  continued  product  times  tho 

analytic  form  from  m  +  1/2  to  *  was  evaluated.  From  tho  table  we  sco 

that  the  analytic  approximation  always  gives  one  percent  accuracy  for 

-20 

reflectivities  larger  than  10 

Evaluation  of  the  continued  product  by  itself  more  often  than  not 
required  thousands  of  factors  to  converge  to  within  10  percent  of  tho 
proper  limit.  Table  A -1  also  presents  some  results  on  this  aspoct .  For 
all  plots  presented  in  Sec.  4  of  this  report,  a  hybrid  calculation  com¬ 
bining  tho  first  10  factors  times  tho  analytic  expression  from  (10  ♦  1  2 


to  infinity  was  used.  We  assume  that  all  points  are  then  accurate  to 
within  one  percent. 

The  form  of  the  dielectric  constant  used  in  Sec.  3  is  a  special 
form  of  that  used  in  Sec.  4.  The  analytic  form  of  Sec.  3  is  obtained 
if  in  Eq .  (A-2)  we  set  o  =  zq  and  let  X  go  to  infinity.  As  a  check  on 
the  validity  of  our  procedures,  we  used  our  analytic  approximation  to 
the  continued  product,  allowed  X  to  approach  infinity  and,  indeed,  found 
convergence  to  the  expression  given  in  Eq .  (6)  of  Sec.  3.  The  procedure 
is  tedious  and  will  not  be  reproduced  here. 
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